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^O Abstract 

r~~i There are many studies dealing with the protection or restoration of wetlands 

\^ and the sustainable economic growth of cities as separate subjects. This study 

\^ investigates the conflict between the two in an area where city growth is threat- 

ening a protected wetland area. We develop a stochastic cellular automaton 
model for urban growth and apply it to the Vecht area surrounding the city 
r^ of Hilversum in the Netherlands, using topographic maps covering the past 150 

'"^ years. We investigate the dependence of the urban growth pattern on the values 

I associated with the protected wetland and other types of landscape surrounding 

^ the city. The conflict between city growth and wetland protection is projected 

On to occur before 2035, assuming full protection of the wetland. Our results also 

^D show that a milder protection policy, allowing some of the wetland to be sacri- 

^^ ficed, could be beneficial for maintaining other valuable landscapes. This insight 

would be difficult to achieve by other analytical means. We conclude that even 
^r slight changes in usage priorities of landscapes can significantly affect the land- 

scape distribution in near future. Our results also point to the importance of a 
protection policy to take the value of surrounding landscapes and the dynamic 
nature of urban areas into account. 
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1. Introduction 

Landscape is defined, in the European Landscape Convention of the |Coun- 



cil of Europe (2000), as "a zone or area as perceived by local people or visitors. 
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whose visual features and character are the result of the action of natural and/or 
cultural (that is, human) factors" . While human activity is one of the actors in 
the formation of landscapes, it is also a cause of their destruction. Today, there 
are unfortunately many landscapes throughout the world that are under seri- 
ous threat due to factors such as intense use of land, pollution, and insufficient 
regional planning. 

Wetlands, the low-lying marshy swampy lands, are a type of landscape under 
direct threat by the current pattern of human development. Even though wet- 
lands have been historically considered as a type of wasteland, their importance 
as integral parts of ecosystems has become increasingly realized during the last 



few decades ( Whigham et al. 1992 ). Now it is known that they strongly support 



wildlife habitats, improve water quality, buffer storms, and control erosion and 



flooding (Kusler and Kentula 1990). Besides their high natural significance, 
they also have a high economic value as a natural resource and as areas for 
recreation. In most developed countries, there are currently many environmen- 
tal management policies in effect to protect wetlands, while trying to permit 
their utilization as a natural resource to continue. 

In this study we intend to analyze a situation where a developing urban 
area lies very close to a wetland under protection. Urban growth is practically 



inevitable in a sustainable economy (Henderson 2005). But can a city resist 



growing to protect a wetland? When will the city threaten the wetland or when 
will the wetland threaten the city? For how long will it be feasible to protect the 
wetland? We investigate and discuss these questions here with a case study of 
the Vecht area in the Netherlands, home to an internationally renowned wetland, 
and the city of Hilversum lying at the border of the protected area (Fig. fll) . 

Our main aim is to investigate the conflicting dynamics of continuing urban 
expansion and the protection of the local wetland. We present urban growth 
projections, based on the current distribution of landscapes and their associated 
importance in terms of natural and economic value, using a stochastic cellular 
automaton (CA) model. The model produces quantitative and visual informa- 
tion on the changing distribution of landscapes in the Vecht area for the near 
future (up to year 2055) and is used to predict the anticipated date of clash 
between the growing city and the wetland. We also study outcomes of several 
different scenarios with different protection priorities. We believe that insights 
from this and similar models can be significantly valuable in short and long term 
regional planning works. 

We provide a brief description of the area in our case study in Section [21 
This is followed by the details of the method, employed computational model, 
and data sources in Section [3] the results of three different runs with the model 
in Section l4J and the conclusions drawn in Section [51 

2. The study area 

2.1. The Vecht area and its urban development policy 

The Vecht area is the floodplain of the river Vecht, forming an important 
part of the delta river system of the Netherlands. It covers a roughly rectangular 
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Figure 1: Overview of the Vcclit area (rectangle), the nearby city of Hilversum, other major 
cities of the Randstad conurbation, and the Green Heart of the Netherlands (shaded area in 
the middle) (Image from Landsat 7. Source: USGS EarthExplorer). 



area of about 8 km by 20 km and includes wetlands, small lakes, streams, and 
fen-grassland patches, which create a characteristic flora and fauna with very 
high biological diversity. 

The history of the area through the last centuries has been a unique combina- 
tion of natural and economic processes: early on dominated by peat extraction 
forming small ponds, transformation of these into larger lakes via erosion, and 
the later recovery of the lost land via drainage and regulation of water tables. 
The area still has a dynamic structure in terms of economy and nature, now 
being trapped in the middle of the Randstad conurbation and considered an im- 
portant part of the National Ecological Network designated by the Netherlands 
Environmental Assessment Agency. 

The Randstad (also known as Deltametropool) is a conurbation in the west- 
ern Netherlands where an increased agglomeration of large cities is underway. 
It is formed mainly by the cities of Amsterdam, Leiden, the Hague, Rotterdam, 
Utrecht, and Hilversum. These cities surround the Green Heart of the Nether- 
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lands, a major open space, of which the Vecht area is a part ( Van Eck et al 
IDG 1997 1 . Fig. fl| gives an overview of the location of the Vecht area rel- 



ative to the Green Heart and the Randstad. While the Green Heart provides a 
valuable natural environment and an attractive leisure area for the population, 
the adjacent cities create new jobs and ask for growing urban development. 
The entire Vecht area is also a part of the "wet axis" of the National Eco- 



logical Network outlined by the Dutch Nature Pohcy Plan (LNV 1997) that 



aims to connect fragmented natural habitats to improve the region's ecology. 



This work, in turn, forms a part of a larger plan, Natura 2000, of the European 
Union for the protection of seriously threatened habitats and the establishment 



of special protection areas throughout Europe ( LNV 2005 ) . 

The simultaneous high demand for urban development and nature preser- 
vation make the task of spatial planning of urban areas in the region highly 
challenging. Although there are planning decisions applied strictly to keep the 
environmental quality high and the development sustainable, progress is yet 
insufhcient as noted in the yearly VINEX reports of the Ministry of Public 



Housing, Spatial Planning and Environment ( |VROM[ [2003|). VINEX states 
that there is a lack of communication and cooperation between municipalities 
(dealing with construction works) and provinces (dealing with the connection 
of natural habitats). 

2.2. The city of Hilversum 

Since Hilversum is located so close to the Vecht area (Fig. [2]), it is immedi- 
ately surrounded by valuable landscapes, comprising not only natural wetlands 
and wetlands arranged as recreation areas, but also the Cornebos Forest, which 



is an old forest attracting rare birds, and the heathland of Westerheide ( VVV 



20061. These areas are highly important in terms of their ecological value. 
Arable lands and pastures nearby also possess economic and ecological impor- 
tance. Pasturelands, to some extent, form a transition zone between the wetland 
and urban areas and provide nesting areas for many bird species. 

Because of the complex structure of the surrounding landscape, city planning 
in Hilversum is strict. The city has managed to maintain a stable population 
for the last 10 years to save the protected landscape from harm, but despite 
low population growth, the extent of the city keeps growing. The number of 
new housings in a region is suggested as an alternative measure of city growth 



by Glaeser and Marshall (1994). Fig. IS] shows the data published by Statistics 



Netherlands (Het Centraal Bureau voor de Statistiek - CBS), indicating a clear 
trend of increase in the number of residential buildings, even though the popu- 
lation remains more or less steady during the same period. This suggests that 
the phenomenon commonly referred to as "urban sprawl" is in effect in the area. 



confirming Kasanko et al. ^2006) 



3. Method 

We develop a stochastic cellular automaton (CA) model for urban growth, 
which is built upon the concept of a total importance valu^ for each landscape 
onto which the city may grow. The model is kept as simple as possible at this 
stage to examine growth dynamics with clarity. We implement the model in 
C# programming language and provide the core source code in the appendix. 
The CA model is further described in Section IX2l 



^A simplistic metric summarizing the economic, natural and other values associated with 
a piece of land. 




Figure 2: The Vecht area and its immediate neighbor, Hilversum (Landsat 7 pseudocolor from 
bands 1-2-4. Source: USGS EarthExplorer). 
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Figure 3: Change in the population and the number of residential buildings in the city of 
Hilversuin during the last decade. Source: Statistics Netherlands - CBS StatLine (available 
online: http://www.cbs.nl/en-GB/menu/cijfers/statline) 



3.1. Input and assumptions 

For making projections with the model, we first need to adopt a reasonable 
average annual rate of urban growth. We assume, by investigating historical 
maps of the Vecht area, one can get a fairly accurate projection of the middle- 
term growth rate, provided that the maps are sufficiently detailed. For this 
purpose, we analyze a collection of topographic maps of the region covering a 
period between 1850 and 1951, with the majority belonging to the 20th century 
(Table IlJ. The reason for not including more recent maps is that during most 
of this selected period there were no restrictions in place to protect natural 
habitats, and as such, the measured growth rate would reflect unrestrained 
factors affecting the growth. This is favorable because it would reproduce the 
unbounded growth pattern of the city, while the limiting pressure of landscape 
protection will be dynamically imposed by the CA model in each simulation 
step on top of this basis growth rate. 

For our measurements, we define the area of the city of Hilversum on each 
historical map by supervised classification using RSI-ENVI 4.1 software. We 
conclude that an annual urban growth rate of around 2% is a reasonable value 
to assume, after observing trends of growth in historical data and giving more 



Table 1: The list of employed spatial data. 



Data 



Source 



Scale / 

resolution 



Date 



Topographic map Dienst der Militaire 1 : 50 000 

of Amsterdam and Verkenningen'^ 

Hilversum 



1850 



Topographic map Topografische Inrichting 1 : 50 000 



1909, 1913, 
1924, 1931 
combined 



Topographic map of Netherlands Ministry of 1 : 50 000 
Utrecht and Amers- War, Topographic Ser- 
foort vice 



1951 



Satellite image GeoCoverTM2000 - 14.25 

Landsat 7 ETM+ m/pixel 
NASA'' 



c. 2000 



Vecht area restora- Van den Bergh et al. 
tion plan 



2001 



2004 



I and Gilbert et al. 
I 



N/A 



2004 



" Provided by Geudeke and Zandvliet j 1990 i^ 

'Accessible online (https : //zulu . ssc. nasa . gov/mrsid ) as compressed MrSID^'^^color 
imagery. 



weight to the trend during the latter half of the 20th century. This value is also 
checked against the short-term projected economic growth rate by Statistics 
Netherlands and is consistent. 

For defining the boundaries of the protected wetland area in our model, we 



employ a suggested restoration plan provided by Van den Bergh et al. ( 2001 ) and 



Gilbert et al. ( 2004 ) , which they use to analyze and evaluate land- use scenarios 



for wetland areas in the Netherlands, but just taking the present situation into 
account. 



3.2. The model 

While being simple discrete systems of recurrently applied rules, cellular au- 
tomaton (CA) models have a demonstrated ability to replicate important fea- 
tures of complexity observed in natural processes, as well as having a number of 
advantages over continuous mathematical models, which are well discussed in 



CA literature (Schiff 2008). Using CA models for describing urban growth has 



been receiving increasing attention in the field of urban planning to study the 



effects of growth and to make predictions in diverse geographic settings ( Itami 



1994 Couclelis 1997 Torrens and O'Sullivan 2001). A well-known CA model of 



this kind is the SLEUTH (Slope, Land cover, Exclusion, Urbanization, Trans- 
portation, and Hillshade) modej^ developed by the United States Geological 
Survey ( USGS 2003[ ) , incorporating factors such as the effect of transportation 



networks and slopes, in addition to the spreading of the existing urban area 



( |Jantz et al.[|2003[ ). 

Here, we believe that a simpler model will be sufficient for the purposes of 
this study. The Vecht area is relatively small compared to the scales urban CA 
models are frequently applied to and this means that the effect of excluding ad- 
vanced features in the model, such as &ree(im(?jj will be insignificant. Another 
important factor is the characteristic flatness of the topography of the Nether- 
lands, where most of the country including the Vecht area is virtually void of 



any hills or slopes (RDG 1979). This renders factors like slope resistance, a 
very important part of general urban CA models, irrelevant in this case. 

In our CA model, we consider the area under study as a grid of cells, where 
each cell is one of 6 types: urban, wetland, heathland, forest, pasture, or water. 
All cell types except the water type can be urbanized. For initializing the model 
with the current distribution of landscapes in the region, we employ a recent 
satellite image (from the GeoCover''''^2000 set distributed by NASA based on 
Landsat 7 ETM-I-, Table IT]) and process it with the supervised classification 
technique. The 6 cell types are assigned real numbers ranging from to 1 
representing their importance, where a value of means no protection at all for 
the cells of that type and a value of 1 means total protection of cells of that 
type under all circumstances. 

Each year, according to the annual growth rate, a number of new cells are 
appended to the existing urban area by converting neighboring non-urban cells 
(i.e. wetland, heathland, forest or pasture) into urban cells. This is done as 
follows: (1) A random cell is selected on the grid, and then checked for whether 
it is non-urban or not, and if it has at least one urban neighbor. The random 
selection continues until such a cell is found. (2) The selected cell is then 
converted into an urban cell, with a probability inversely proportional to the 
landscape importance value assigned to its type. (3) For every year during 
the simulation, this random selection and conversion cycle continues until the 
required number of new urban cells for that year are created. The model permits 
the formation of diverse development scenarios by assigning different importance 
values to different landscape types. The flowchart in Fig. |4] gives an outline of 
the model's working. 

The number of new urban cells created each year is simply 

N' = rx N (1) 

where N' is the number of new urban cells, r is the annual growth rate, and 
N is the number of present urban cells. The annual growth rate of urban area 



■^http : //www . ncgla . ucsb . edu/proj ects/gig/' 



''a type of growth in urban CA models, in which a new isolated urban cell is occasionally 
formed without any adjacent urban cells 



is held fixed during a run, based on the assumption that the expected growth 



rate of a city should be independent of its size (Cordoba 2004). 



Partition the map into urban, wetland, heatliland, 

forest, pasture and water cells, 

set on a rectangular grid. 



Assign relative landscape importance values to 
different cell types. 
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Assume an urban growth rate. 



Calculate the number of new urban points N, 
according to growth rate. 
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with a probability inversely proportional to its 

importance value. 
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Repeat until A^ non-urban cells are converted. 
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Repeat these for the next year, until the final year. 



Figure 4: Flowchart of the employed urban growth model. 



We take a cautious approach by presuming an ultimate date of reasonable 
prediction for this model before the year 2055, thinking that calculations cov- 
ering longer periods would become more diverted by the accumulation of errors 
introduced by assumptions, selected parameters such as the growth rate, and 



the intrinsic limitations of the CA modeling approach (Yeh and Li 2006). Also 



the trends and priorities in the Vecht area and in the Netherlands in general are 
likely to be changed beyond terms longer than this. 



4. Results and discussion 

Going back to the questions at the beginning: 



4--1- Can the city resist growing to protect the wetland? 

This part cannot be answered solely with the aid of the model, but literature 
and experiences may help. Economic growth is always desired and has a direct 



influence on the growth of cities. A similar situation described by Karen (2001 1 
in Merced County in the United States, an area that has to deal with urban 
growth and wetland recreation simultaneously in California's Central Valley, 
exemplifies how complex the situation could become. It has been noted that $ 
27.7 million per year is spent for the maintenance of the wetland recreation area, 
and the wetland contributes $ 41 million per year back to the local economy and 
provides jobs for more than 800 people. The paradox posed by this situation 
can immediately be seen: the effort to protect the wetland turns back as a 
contribution to the local economy, which in turn promotes urban growth and 
threatens the wetland. 

At the same time, wetlands are now recognized globally as a cornerstone 
and focal point of economic development in both the developed and developing 
countries and ways are sought after that can balance their protection and use. 
Wetland restoration in the Vecht area presents many valuable economic benefits 
to nearby cities, such as income from increased tourism and fishing; and for 
the city of Hilversum, economic growth is foreseeable. This growth should 
be planned in a way which ensures the survival of the wetland and the other 
landscapes in the region. 

4-2. When will the city threaten the wetland or when will the wetland threaten 
the city? 

To answer this question, we design the first run with the model using param- 
eters thought to represent current policies / priorities in the region. Wetland 
importance value was set to 1.00 (i.e. it is impossible to destroy wetland) to 
make sure that it will be protected under all circumstances. Used parameter 
values are summarized in Table [5] 

The main result from our simulation, in addition to the projected spatial 
distribution on the map, is the amount of loss in the area of each landscape type 
present. Here, the resulting distribution of different landscapes is summarized in 
Table [2] for years 2030 and 2055 (as the amount of area lost and as a ratio of the 
projected area to the current area) and the projected distribution of landscapes 
are presented in Fig. [5] from 2015 to 2055 with 20 years interval. Note that the 
annual urban growth rate during all runs is assumed as 2.0 %, as explained in 
Section |3l 

We predict that the clash between the wetland area and the city bound- 
ary will become significant by the year 2035, well within a timescale to grant 
considering in the current recreation and urban plans of the Vecht area. 

4.3. How long is this wetland going to be protected? 

The results of the first run in the previous section indicate a deficiency with 
that scenario. While the wetland is successfully protected, this comes at the 
cost of losing other valuable landscapes. The most noticeable is the decrease in 
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Figure 5: The first run of tlie model with definite protection for the wetland and decreasing 
importance values for heathland, forest and pasture (details in Tablelll. Resulting distribution 
of landscapes is shown between 2015 and 2055 with 20 year intervals. 



Table 2: Parameters and results for the first run. Definite protection for wetland, decreasing 



importance for heathland, forest and pasture". 



Landscape Importance 2030 2055 

Amount Remaining Amount Remaining 

lost (km^) (%oforig- lost (km^) (% of orig- 
inal) inal) 



Wetland 


1.00 





100.0 





100.0 


Heathland 


0.95 


1.39 


74.58 


3.62 


34.03 


Forest 


0.80 


5.01 


88.22 


14.47 


66.01 


Pasture 


0.65 


10.17 


90.82 


26.64 


72.95 



" Annual growth rate is 2.0 %. The simulation starts from year 2005. 

the total forest area to almost 66% of its original amount by 2055 (Table [2]). 
One may argue that the proximity of the forest to the urban area makes it more 
prone to depletion, but it should still be possible to prevent an amount of this 
loss by trying to adjust the protection priorities. 

The other two runs of the model in this section are designed to investigate 
whether permitting some loss of wetland helps to save a greater amount of 
the forest. With this purpose, the model was run many times, giving different 
importance values to each type of landscape and observing the results. Here we 
present results from two such runs. 

Table|3]and Fig. [6] present a situation where the wetland importance param- 
eter is reduced to 0.90, in an effort to prevent the drastic loss of forest in the 
first run. Surprisingly, it is noticed that with a little sacrifice of the wetland, one 
has the potential to save a significant amount of forest: around 82% remaining 
by 2055, in contrast to the figure of 66% in the first run (compare Table [2] and 
Table [3|. 

The results of a more radical change in protection policy are presented in Ta- 
ble |4] and Fig.JTJ where we assign wetland, heathland and forest equal protection 
priorities. This results in a situation in which a more even distribution of all 
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Figure 6: The second run of the model, with the wetland importance reduced to 0.90, in an 
effort to prevent the drastic loss of forest in the first run (details in Table pi. The resulting 
distribution of landscapes is shown between 2015 and 2055 with 20 year intervals. 



Table 3: Parameters and results for the second run. The second run with the model. 



Landscape 


Importance 




2030 


20 


55 






Amount 




Remaining 


Amount 


Remaining 






lost (km' 





(% of orig- 
inal) 


lost (km^) 


(% of orig- 
inal) 


Wetland 


0.90 


0.87 




97.30 


2.46 


92.39 


Heathland 


0.90 


1.69 




69.27 


4.35 


20.68 


Forest 


0.85 


2.75 




93.53 


7.69 


81.93 


Pasture 


0.65 


9.13 




91.76 


23.55 


78.75 



" Annual growth rate is 2.0 %. The simulation starts from year 2005. 



resources is achieved while still maintaining the wetland above others. An issue 
to recognize here is how this is achieved while the wetland is assigned the same 
importance value as the other landscapes. This is the result of the particular 
geographical distribution of landscapes in this region: forest and heathland are 
under immediate threat by the growth of Hilversum, while the wetland lies at 
some distance, as compared to these. 

If wetland protection is considered of highest importance in the Vecht area, 
then the city of Hilversum will totally grow on other landscapes. In this case, 
these landscapes will be in danger of irreversible loss to keep the wetland safe. 
The results presented here indicate that setting a conventional, absolute con- 
straint with the purpose of protecting a type of landscape can lead to other 
serious problems as a side effect. We conclude that an environmental man- 
agement project should take the complex interaction and the unique spatial 
arrangement of different landscapes in a region into account, instead of sim- 
ply putting particular regions into protected status according to the immediate 
perceived threat of the day. 
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Figure 7: The third run of the model with equal importance for the wetland, heathland and 
forest (details in Table |4]|. Resulting distribution of landscapes is shown between 2015 and 
2055 with 20 year intervals. 



Table 4: Parameters and results for the third run. Equal importance for wetland, heathland 
and forest". 

Landscape Importance 2030 2055 

Amount Remaining Amount Remaining 

lost (km^) (%oforig- lost (km^) (% of orig- 
inal) inal) 



Wetland 


0.80 


1.16 


96.40 


3.38 


89.54 


Heathland 


0.80 


2.28 


58.36 


4.35 


20.65 


Forest 


0.80 


2.33 


94.52 


6.95 


83.67 


Pasture 


0.65 


8.69 


92.15 


23.29 


78.98 



" Annual growth rate is 2.0 %. The simulation starts from year 2005. 

4.4- Further discussion 

There is another way to protect those landscapes from the threatening urban 
growth in Hilversum, and it is even more effective: promoting a vertical growth 
(i.e. higher buildings) instead of the current trend of lateral growth. This 
solution can be a very efficient in protecting the valuable landscapes in the 
Vecht area, but there is no indication that this will be happening in the near 
future. 

An important insight gained from this study is on the mode of propagation 
of a growing city in such a configuration of landscapes and the speed at which 
this occurs. This type of information can be very useful in environmental man- 
agement projects by making it easier to pick which urban areas in the proximity 
of a protected area to get included in the protection plan (by predicting whether 
an urban area will pose a significant threat during the course of the protection 
plan). 

Another issue of importance is how the suggested protection priorities, as 
represented in our model by the so-called importance values, could be imple- 
mented in reality, while how these priorities should be determined remains a 
separate subject. The issue often exceeds environmental considerations and in- 
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volves other factors such as community demand for each landscape type. A 
recent study by |Kaplan and Austin ( 2004 ) on community preferences about 



neighboring different landscapes, for instance, has shown that there is a consid- 
erably greater desire for forests than other landscapes, regardless any benefits 
that might be obtained from these. 

5. Conclusions 

Although urban growth is a dynamic and complex process, significant in- 
sight can still be gained by experimentation with computational models. An 
investigation of the Vecht area and the city of Hilversum is particularly useful 
in demonstrating this modeling and prediction approach, given the diversity of 
landscapes present in the region and the immediate threat posed on these by 
urban growth. 

We predict that the pressure to protect the designated wetland area will 
probably keep the city away from the wetland until around 2035, after which 
the expected clash between the city and the protection area will occur. But up to 
that time, the city would have to grow into other valuable landscapes, including 
forests. This leads to our suggestion that, to protect not only the wetland but 
also the diversity of other landscapes, the protection pressure on the wetland 
should be reduced to a level bringing the anticipated clash to an earlier date, 
while still preserving a considerable amount of wetland. This way, we believe, 
a more even distribution of all resources is achieved while still maintaining the 
wetland above others. An important aspect to recognize here is how this can be 
achieved even while the wetland is assigned the same importance value as the 
other landscapes. 

It is certain that landscapes will change continuously. But this change should 
be monitored and controlled in a way that is most beneficial for nature and hu- 
man society. By investigating this interplay through experiments with different 
protection scenarios, this study shows that even a slight adjustment in the uti- 
lization priorities of landscapes is capable of achieving important changes in the 
distribution of landscapes that we will be depending on in the near future. 
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Appendix A. Code for the CA model 

Following is the core C# code of the CA urban growth model used in this 
study. For brevity, only the portion central to the working of the model is 
presented, leaving out auxiliary code such as for the interface. The code was 
compiled and executed on Microsoft .NET Framework version 2.0.50727. 

namespace HilversumVecht 

{ 

public class CAModel 

{ 

private enum CellType {Forest, Heathland , Urban, 
Wetland, Water, Canal, Void} 

//Predefined RGB values for: 

//I. Reading the cell types from the initial map 

presented to the class 
//as a bitmap via the constructor (Bitmap inttialMap). 
//2. Producing the output bitmap while the model is 

running . 
private static Color ModelMapForestCoIor = Color. 

FromArgb ( 1 2 1 , 174, 4); 
private static Color ModelMapHeathlandCoIor = Color. 

FromArgb ( 2 1 9 , 198, 215); 
private static Color ModelMapUrbanAreaColor = Color. 

FromArgb(231, 56, 13); 
private static Color ModelMapWetlandColor = Color. 

FromArgb(13, 111, 55); 
private static Color ModelMapWaterColor = Color. 

FromArgb (46, 143, 210); 
private static Color ModelMapCanalColor = Color. 

FromArgb(107, 175, 220); 
private static Color ModelMapVoidColor = Color. 

FromArgb(195, 212, 144); 
private CellType [,] ModelMap ; 
private int Width, Height; 
private float ForestValue , WetlandValue , 

HeathlandValue , PastureValue ; 
private bool Abort; 
private Random rnd ; 

// The constructor of the class . 

public CAModel( int width, int height , Bitmap 

initialMap , float forestValue , float wetlandValue , 
float heathlandValue , float pastureValue ) 

{ 

Width = width ; Height = height ; ForestValue = 
forestValue ; WetlandValue = wetlandValue ; 
HeathlandValue = heathlandValue ; PastureValue 
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= pastureValue ; 
rnd — new Random () ; 
ModelMap = ReadModelMap( initialMap ) ; 

} 

// Reads the initial landscape distribution from the 

bitmap . 
private CellType [ ,] ReadModelMap( Bitmap map) 

{ 

CellType[,] ret = new CellType [Width , Height]; 
for (int X = ; x < Width ; x++) 

for (int y = ; y < Height ; y++) 

rct[x, y] = ColorToCellType(map. GetPixel(x 

, y)); 

return ret ; 
} 

// Converts a color in the initial landscape 
distribution bitmap to a CellType value. 
private static CellType ColorToCellType ( Color c) 

{ 

if (c ^= ModelMapForestColor) 

return CellType . Forest ; 
else if (c ^= ModelMapHeathlandColor ) 

return CellType . Heathland ; 
else if (c = ModelMapWetlandColor) 

return CellType . Wetland ; 
else if (c = ModelMapUrbanAreaColor) 

return CellType . Urban ; 
else 

return CellType . Void ; 
} 

// Converts a CellType value into the corresponding 

color for producing the bitmap output. 
private static Color CellTypeToColor ( CellType c) 

{ 

if (c ^= CellType. Forest) 

return ModelMapForestColor ; 
else if (c = CellType . Heathland ) 

return ModelMapHeathlandColor ; 
else if (c ^= CellType . Wetland) 

return ModelMapWetlandColor; 
else if (c = CellType . Urban) 

return ModelMapUrbanAreaColor ; 
else if (c = CellType . Void) 

return ModelMap VoidColor ; 
else 

return Color . Yellow ; 
} 
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// Main method for running the model. 

public void Run( float growthRate , int steps) 

{ 

Abort = false ; 

float init ialForest Area = CountArea( CellType . 

Forest ) ; 
float initialWetlandArea = CountArea( CellType . 

Wetland) ; 
float initialHeathArea = CountArea( CellType . 

Heathland ) ; 
float initialUrbanArea = CountArea( CellType . Urban) 



} 



float initialOther Area = Width * Height — 

initialUrbanArea — initialWetlandArea — 

initialForest Area ; 
int fa , wa, ha , ua ; 
float percentForest Area , percentWetlandArea , 

percentUrbanArea , percentHeathArea , 

percentOtherArea ; 
Bitmap map = new Bitmap (Width , Height) ; 
for (int i = ; (i < steps) && ! Abort ; i++) 

{ 

UpdateModelMap( growthRate) ; 
fa = CountArea( CellType . Forest ) ; wa = 
CountArea( CellType .Wetland) ; ha = 
CountArea( CellType . Heathland) ; ua = 
CountArea ( CellType . Urban) ; 
percentForest Area = (fa / init ialForost Area ) * 
100; percentWetlandArea = (wa / 
initialWetlandArea) * 100; 
percentUrbanArea = (ua / initialUrbanArea) 

* 100; percentHeathArea = (ha / 
initialHeathArea) * 100; percentOtherArea 
= ((Width * Height — fa — wa — ua — ha) / 
initialOther Area ) * 100; 
map = WriteModelMapO ; 
} 



// The annual update to the grid. 

private void UpdateModelMap( float growthRate) 

{ 

int px , py , timeout ; 

int steps = ( int )(( growthRate / 100) * CountArea( 

CellType. Urban)) ; 
for (int i = ; i < steps ; i++) 

{ 

//Timeout introduced for performance concerns 
timeout = 100000; 
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while (timeout > 0) 

{ 

px = rnd.Next(l, Width); py = rnd . Next ( 1 , 

Height); 
if (ModelMap[px, py] != CellType . Urban) 

{ 

if ( UrbanNcighbors (px , py) > 2) 

{ 

if ( ModelMap [ px , py] = CellType. 
Void) 

{ 

if ( PastureValue <= rnd. 
NextDoubleO ) 

{ 

ModelMap [ px , py ] = 

CellType . Urban ; break; 
} 
} 

else if (ModelMap [px , py] ^= 
CellType. Forest) 

{ 

if (ForestValuc <= rnd. 
NextDoubleO ) 

{ 

ModelMap [ px , py ] = 

CellType . Urban ; break; 
} 
} 

else if ( ModelMap [ px , py] ^= 
CellType. Wetland) 



{ 



} 
else 

{ 



if ( WetlandValue <= rnd. 
NextDoubleO ) 

{ 

ModelMap [ px , py ] = 

CellType . Urban ; break; 
} 



} 
} 
} 
timeout 



if ( HeathlandValue <= rnd. 
NextDoubleO ) 

{ 

ModelMap [ px , py ] = 

CellType . Urban ; break; 
} 
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} 

// Produces the output bitmap from the current grid ( 

variable "ModelMap" ) . 
private Bitmap WriteModclMap () 

{ 

Bitmap ret = new Bitmap( Width , Height); 
for (int X = ; x < Width ; x++) 

for (int y = ; y < Height ; y++) 

ret . SetPixel(x, y, IntegerToCoIor (ModelMap [ 

X, y])); 

return ret ; 

} 

// Counts the number of urban neighbors of a given 

cell . 
private int UrbanNeighbors ( int x, int y) 

{ 

//Moore neighborhood 
int ret = 0; 



if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
if (ModelMap 

ret++; 
return ret ; 



X - 1, y - 1] ^ CellType. Urban) 

x, y - 1] ^ CellType. Urban) 

X + 1, y - 1] ^ CellType. Urban) 

X - 1, y] ^= CellType .Urban) 

X + 1, y] ^ CellType .Urban) 

X - 1, y + 1] ^ CellType. Urban) 

X, y + 1] ^ CellType .Urban) 

X + 1, y + 1] ^ CellType .Urban) 



} 



// Counts the total number of cells of a given 

CellType on the grid. 
private int CountArea( CellType c) 

{ 

int ret = 0; 

for (int X = ; x < Width ; x++) 

for (int y = ; y < Height ; y++) 
if (ModelMap [x, y] ^ c) 
ret++; 
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return ret ; 
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